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Abstract 

Asymptotic expansions are given for large values of n of the gener- 
alized Bernoulli polynomials B£(z) and Euler polynomials E£(z). In 
a previous paper Lopez and Temme (1999) these polynomials have 
been considered for large values of fi, with n fixed. In the literature 
no complete description of the large n asymptotics of the considered 
polynomials is available. We give the general expansions, summarize 
known results of special cases and give more details about these re- 
sults. We use two-point Taylor expansions for obtaining new type of 
expansions. The analysis is based on contour integrals that follow from 
the generating functions of the polynomials. 
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1 Introduction 

Generalized Bernoulli and Euler polynomials of degree n, complex order \x 
and complex argument z, denoted respectively by Bn(z) and En(z), can be 
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defined by their generating functions. We have [5, 10] 



-W , \W\ < Z7T, (1.1) 



v ; n=0 

and 

_Zf^_ = f «M tt ,«, w< ,. (1 . 2) 

(e w + iy ^ n! 1 1 1 ; 

v ' n=0 

These polynomials play an important role in the calculus of finite dif- 
ferences. In fact, the coefficients in all the usual central-difference formulae 
for interpolation, numerical differentiation and integration, and differences 
in terms of derivatives can be expressed in terms of these polynomials (see 
[5,7]). 

An explicit formula for the generalized Bernoulli polynomials can be 
found in [9]. Properties and explicit formulas for the generalized Bernoulli 
and Euler numbers can be found in [4, 11, 12] and in cited references. 

In a previous paper [1] we have considered these polynomials for large 
values of fi, with n fixed, and in the present paper we consider n as the large 
parameter, with the other parameters fixed. We summarize known results 
from the literature for integer values of (i, and give more details about these 
results. We describe the method for obtaining the coefficients in the expan- 
sion for general \i. Finally, we use two-point Taylor expansions for obtaining 
new type of expansions for general [i. The analysis is based on contour in- 
tegrals that follow from the generating functions of the polynomials. 

2 The generalized Bernoulli polynomials 

Three different cases arise, depending on jjl = 0, —1, —2, . . ., [i = 1, 2, 3, ... , 
and n otherwise, real or complex. Our approach is based on the Cauchy 
integral 

ti\ f w^& wz dw 

where C is a circle around the origin, with radius less than 2ir. This follows 
from (1.1). 

2.1 Asymptotic form when fi = 0, —1, —2, . . . 

In this case the generating series in (1.1) converges for all finite values of 
z, and, hence, the polynomials Bn(z) have a completely different behavior 
compared with the general case. 
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We first follow the approach given in [14], and observe that when \i is 
a negative integer or zero, say fx = —m (m = 0,1,2,...), we can express 
Bn(z) in terms of a finite sum. We expand by the binomial theorem 



TO , <. 

(e w - l) m = ^(-i) m "M m ), 



(2.2) 



This gives 



(n + m) 



m , 
r=0 ^ 



(z + r) 



n+m 



(2.3) 



For any given m £ N and complex z only the term or terms with the largest 
values of \z + r\ will give a large contribution to the sum in (2.3), the other 
terms being exponentially small in comparison. We conclude, that (2.3) 
gives the asymptotic form when n — > oo, when fi = —m and z are fixed. 

In particular, when z > 0, the term with index r = m is maximal, and 
we have 



B- m (z) 



(n + m)! 



(z + m 



n+m 



l + O 



. -i \ n+m 

Z + m — 1 \ 

z + m J 



(2.4) 



The error term can also be estimated by C(exp(— (n + m)/(z + m))), which 
is indeed exponentially small compared with unity. 

For general complex z = x + iy and x > —m/2 the term with index 
r = m again is maximal and the same estimate as in (2.4) is valid. When 
x = —m/2 the terms with r = and r = m give the maximal contributions, 
and we have 



B~ m {z) 



(-l) m (-im + iy) n+m + am + zy) 



n+rn 



(n + m)! 

When x < —m/2 the term with index r = is maximal, and we have 



(2.5) 



B~ m (z) = (-iy 



nl 



n+m 



(n + m)! 



l + O 



z + l 



n+m 



(2.6) 



By using the saddle point method we can obtain an estimate similar as 
the one in (2.4). We write 



B r 



nl 



[1 - e-) 



dw 



w 



(2.7) 
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where 4>{w) = (z + m)w — (n + m) hiw. This function has a saddle point 
at wo = (n + m) / (z + m), and when $two is large and positive we replace 
(1 — e~ w ) by its value at this point. Then we have 



B- m (z) ~ h- e - w °) m — [ e^ m > 

" 1 ^ V ' 2-Ki Jc w n+m+l 

v 7 (n + m)! v 7 



(2.8) 



in which e _u, ° is exponentially small. When 3ftu>o is not positive we can 
modify this method to obtain the estimates as given in (2.5) and (2.6). Also, 
we can make further steps in the saddle point analysis, and show that the 
next terms in the expansion are exponentially small compared with unity, 
similar as shown in (2.4). However, the representation in (2.3) describes 
very elegantly the asymptotic behavior. 

2.2 Asymptotic form when /n = 1, 2, 3, . . . 

We write fi = m. The starting point is the expansion for m = 1: 

oo 2'Kikz 

Hiw = b„(,) = -«! E £sir . w 

k — — oo v 7 
fe^O 

which for z G (0, 1) can be viewed as a Fourier expansion of B^(z). This 
expansion follows from taking the radius of the circle C in (2.1) equal to 
(2K + l)ir, {K an integer). Taking K large, we take into account the poles 
of the integrand at w = 2-nik (k = ±1, ±2, . . .), and calculate the residues 
of these poles. The integral around the circle C tends to zero as K — > oo, 
provided n > 1 and < z < 1. 

This gives the expansions (n = 1, 2, 3, . . . ; < z < 1) 

ft.W-2(-l)-«(».)lf S^^, (2.10) 

and 

B aB+1 (z) = 2(-irH2n+l)!E^PS- 

In (2.11) we can take n = 0, provided < z < 1. This gives the well-known 
Fourier expansion of B\(z) = z — 1/2, < z < 1. 
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In (2.9) only the terms with k = ±1 are relevant for the asymptotic 
behavior, and we obtain for fixed complex z 

B l n {z) = { J • [cos (27TZ + \im) + O (2~")j , n ^ oo. (2.12) 

For general fixed real or complex z the series in (2.10) and (2.11) can be 
viewed as asymptotic expansion for large n, as easily follows from the ratio 
test. 

For general \i = m = 1, 2, 3, . . . a similar expansion as in (2.9) can be 
given. In that case the poles of the integrand in (2.1) are of higher order. 
We can write 

oo 2nikz 

B™(z) = -n\ £ PT(n,z)^—-, (2.13) 

fc— — oo * ' 

fe^O 

where (3l(n,z) = l,Vfc. This is a Fourier expansion for z £ (0,1) when 
m < n. For other values of z it can be used as an asymptotic expansion for 
large n. 

An explicit form of /3™(n, z) follows from calculating the residues of the 
poles at 2irik of order m of the integrand in (2.1). For this we compute the 
coefficient c m _i in the expansion 

(w - 2nik) m w m e zw ^ , n , ft1i , 
— — = > cJw - 2mk) r , (2.14) 

V 7 r=0 

from which we obtain 

^hw' =c -- 1 - (2 - 15) 

We substitute w = s + 2-Kik and write the expansion as 

= ^c r s r . (2.16) 



e 



a m e zs 
2-nikz b e 



(e s - l) m (s + 2vrifc) n + 1 - r 



We use (1.1) and write the left-hand side in the form 



(2mk) n+1 ~ m ^ n\ ^ V i/ 7 (2vri/c)^ 
Hence, c r of (2.16) can be written as 



e 



(2*i*)' 
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and we conclude that 



Jl-wikz ™ 1 Dm/ 



Cm-1 



B™(z) fm — n — 1 



(2iTik) n ^— ' i/! \m — 1 — v 



(2irik) v . 



It follows from (2.15) that 



m— 1 



^=0 



B™(z) fm — n — 1 
m — 1 — i/ 



{2-Kikf. 



(2.19) 



(2.20) 



To avoid binomials with negative integers, and to extract the main asymp- 
totic factor, we write 



m—l 



n — 1 

m 



m—l 



m — 1\ (n — ^ — 1)! 



v ) (n-1)! 



For large n the main term occurs for v = 0. We have 

(_ 1 jm-l n m-l 



/3T(n,z) = 



[l + OCn" 1 )]. 



-2mk) v . 
(2.21) 

(2.22) 



(m-1)! 

Observing that, as in (2.9), only the terms with k = ±1 are relevant for 
the asymptotic behavior, we obtain 

(_l)n+l„! 



(2vr) n 



/?r(n, z)e 2 ^+ 2™ + #™ (n, *)e 



2"*'° + . . . , 

(2.23) 



and by using (2.21) we obtain for fixed m and complex z (cf. (2.12)) 
2(_!)m+n f n -r 



m—l 



(27T) 



i/=0 



m — l 
m — 1\ (n — i/ — 1)! 



i/ J (n-1)! 



(2vrf cos a + O {2- n ) 



(2.24) 



as n — > oo, where <r = (2z + — j^)vr. 

To obtain /3™(n, z) for m > 1 we can also use a recurrence relation. We 
have the relation 



HB£ +I (z) = (u - n)Btt(z) + n(z - rfB^z), n > 1, 



(2.25) 



which follows from (1.1) by differentiating both members with respect to w. 
By differentiation with respect to z we find 



nB^iz) = j- z B»{z 



(2.26) 
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Figure 1: Contour for (2.1) for general \i 



371/ 




giving 

f,B%+\z) = (fi-n)B£(z) + (z-»)-?-B£(z), n > 0. (2.27) 
This gives the recurrence relation for m = 1, 2, 3, . . . 
m(3™ +1 (n,z) = [m-n + 2mk(z-m)]/3^(n,z) + (z-m)-^f3^(n,z). (2.28) 

2.3 Asymptotic form for general complex \x 

We consider (2.1) and observe that the singularities at ±2ni are the sources 
for the main asymptotic contributions. We integrate around a circle with 
radius 3tt, avoiding branch cuts running from ±2ni to +oo. See Figure 1. 
The contribution from the circular arc is 0((37r)~ n ), which is exponentially 
small with respect to the main contributions. 

We denote the loops by C± and the contributions from the loops by I±. 
For the upper loop we substitute w = 2nie s . This gives 



where 

/ 27T? s \ M 

^ = [-^—j) eZU+>iS > u = 27ri(e s -l), (2.30) 

and C_|_ is the image of £+. C + is a contour that encircles the origin in the 
clockwise fashion. 

To obtain an asymptotic expansion we apply Watson's lemma for loop 
integrals, see [8, p. 120]. We expand 

oo 

g(s) = Y J 9kS k , (2.31) 

k=0 

substitute this in (2.29), and interchange summation and integration. This 
gives 

e 2wiz 00 

where 

Fk = 1 tdl sk ^ e ~ nSds > ( 2 - 33 ) 

with C + extended to +oo. That is, we start the integration along the contour 
C + at s = +oo, with phs = 2ir, turn around the origin in the clock- wise 
direction, and return to +oo with phs = 0. 

To evaluate the integrals we turn the path by writing s = e m t, and use 
the representation of the reciprocal gamma function in terms of the Hankel 
contour; see [10, p. 48]. The result is 

k ~ n 6 r(/x - k) ~ r( M ) ' 1 j 

where (a) n is the shifted factorial, or Pochhammer's symbol, defined by 

T(a + n) 

(a) n = a ■ (a + 1) • • • (a + n - 1) = V ' , n = 0,l,2,.... (2.35) 

1 (^aj 

This gives the expansion 

nln^ 1 , ^ (1 - ^) fc g fc 
(2vr)«r(^) ^ n fe 

where 

X = 2C - V, C = (^ + h)7r. (2.37) 
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The contribution /_ can be obtained in a similar way. However, it is 
the complex conjugate of I + (not considering z and fi as complex numbers). 
The result i+ + 7_ can be obtained by taking twice the real part of I+. We 

(r) (i) / (r) (i) 

write gk = g k + ig k (with g k ,g k real when 2 and fx are real) , and obtain 



2 77' r?^" 1 



(27r)»r( A1 ) 



cosx V e smx > e 

Z — / n K L — / n « 

k=0 k=0 



(2.38) 



as n — > 00, with 2 and /U fixed complex numbers (// ^ Z). 
The first few coefficients gjf , g^ k are 

(r) 1 (») n 

^ = i/i, <?? = 2C, 

= U^ 2 + (4vr 2 - IV - 48C 2 ), = (1 + M)C, 

4° = A(M 3 + (4vr 2 - 1)^ + 8(^ 2 - QQ^ - 96C 2 ), 
= AC(3^ 2 + (4vr 2 + 5V-16C 2 + 4). 



(2.39) 



The first-order approximation reads 



2n'n^ _1 r 1 

mz) = Wrnd [ c ™< 2z +v- 2 n ) + °( l i n )\ > ( 2 - 4 °) 

Norlund [6, p. 39] describes the same method of this section and only gives 
the first-order approximation. 

2.3.1 An alternative expansion 

As observed in the previous method, the main contributions to (2.1) comes 
from the singular points of the integrand at ±27ri. In this section we expand 
part of the integrand of (2.1) in a two-point Taylor expansion. In this way 
a simpler asymptotic representation can be obtained. For more details on 
this topic we refer to [2, 3] and for the evaluation of coefficients of such 
expansions to [13]. We write 

f(w) = {w 2 + 4vr 2 ) M (^y) " e wz (2.41) 
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and expand 

oo 

f(w) = K + v>0 k ) {w 2 + 47T 2 ) k . (2.42) 

fc=0 

The function f(w) is analytic inside the disk \w\ < 4tt and the series con- 
verges in the same domain. The coefficients Qo and /3q can be found by 
substituting w = ±2iri. This gives 

f(2m) + f(-2m) 
q = = cos2C, 

_ /(27TQ - f(-2m) _ 1 ^ 
A_ 4™ -2^ Sm2C ' 

where ( is defined in (2.37). 

The next coefficients can be obtained by writing fo(w) = f(w) and 

/jM-(aj + w/3j) 

^ +lH " ^Ti^ 

^ t • i 2 - 44 

= £ (a fc + u,&)(«, 2 + 47r 2 ) fc - J -\ 

fc=i+i 

j = 0, 1, 2, . . ., and by taking the limits when iu — > ±27ri We have 

/j(2™)-/j(-2™) 



&+1 = " 



JJ + ' ~ 167T 2 

This gives 
ai = - — ^ [3/U cos 2( + 2m] sin 2(] , 

lD7T z 

/?i = ^ [2tt»j cos 2C + (2 - 3/i) sin 2(] , 

a 2 = 15367r 4 [(^12tt 2 7? 2 + 4^tt 2 - 33^ + 27^ 2 ) cos 2( + 12ttt/(3^ - 1) sin 2C] , 

02 = e [-36777/^ - 1) cos 2( + (36 - 69/i + 27,u 2 + 4/wr 2 - 12tt V) sin 2(] , 

30727r D 

(2.46) 

where rj = \x — 2z. 

Substituting the expansion in (2.42) into (2.1) we obtain 

oo 

B£{z) = n\ 2 V £ [«*** B) + A** B " 1) ] > ( 2 - 47 ) 

k=0 
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where 



fc 2tu 7 C v ; w n+1 v 7 

We have $j; 2n+1) = and 

$( 2 ™) = ( 27r )2*-2M-2n ^ - ^ = (_l)»(2j)M-2^Ml. (2.49) 



Hence, 

oo 

<(,) = (2n)!2 3 ^ 2 ^a fc $f), 



fc=o 

oo 



(2.50) 



i?L + iW = (2n + l)!2 3 ^^/3 fc <I>f). 

fc=0 

These convergent expansions have an asymptotic character for large n. This 
follows from (see (2.35)) 



^ (/i " k) n 

A; 

^ fc- l + n) r(/i-fc) (2.5i) 

r(/i-ife-i) r(/i-ife + n) 

= 4vr 2 = 0( n - 1 ), n-oo. 

/i — k + n — 1 

We compare the first term approximations given in (2.40) and those from 
(2.50). From (2.40) we obtain 

BL(z) ~ (-1) {2 *yn m cos < 2z + M) + • • • , (2-52) 
and from (2.50) 

Because r(n + /x)/n! ~ n M_1 asm oo, we see that the first approximations 
give the same asymptotic estimates, and they are exactly the same when 
Ai = l. 
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Integer values of /i, Comparing the expansions in (2.38) and (2.50), we 
observe that those in (2.50) do not vanish when ji = 0, —1, —2, . . ., whereas 
the expansion in (2.38) does. We have when fi = m (integer) 



k 



( (27r) 2fc-2m-2n ^ k >n + m 



(2.54) 

0, k < n + m. 



Hence, the summation in (2.50) starts with k = n + m. The scale {&f n } 
loses its asymptotic property, because now 

^) = 4?r ^ + i = H > ( 2 - 55 ) 

where /c = n + m + £, and a possible asymptotic character of the series 
in (2.50) has to be furnished by the coefficients otkifik-, which depend on n 
when k > n + m 

Because n is assumed to be large, and the coefficients a^^k m (2.50) 
become quite complicated when k > n + m, these expansions are of no use 
when fi is an integer. 

When we replace the expansion in (2.42) with 



f(w) = £ (5* + wK) W ^ (2-56) 
fc=o \ w J 



we obtain 



r(2n) 

fc=0 

< +1 (z)~(2n + l)!2 3 ^^$ 



(2n) 
fc 

fe=0 



(2-57) 



where and can be obtained from a similar scheme as in (2.44)-(2.45). 
The functions 3*)?"^ are given by 



-(2n) = (-1)"+* (tl-k) n+k = (-lf+fc r(/* + n) 

fe ( 2vr )2 M +2n ( n + jfe)i (27r) 2 ^+ 2n r(//-fc)(n + fc)!' 



When /j, = m (integer) these functions vanish if k — m = 0, 1, 2, . . ., which is 
more useful that in the earlier choice (2.42). When m < all terms vanish, 
when m > the series have a finite number of terms. 
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With the expansion in (2.56), which converges in certain neighborhoods 
of the points w = ±2ni (and not in a domain that contains any allowed 
deformation of the curve C in (2.1)), the expansions in (2.57) do not converge, 
but they have an asymptotic character for large n. 

As an example, when m = 1 the expansions in (2.57) have just one term 
(k = 0). In this case ®f n) = (-l) ri /(2vr) 2n+2 and 5 = a , A) = A) (see 
(2.43)). These approximations correspond exactly to the first terms in the 
expansions in (2.10) and (2.11). 



3 The generalized Euler polynomials 

We can use the same methods as for the Bernoulli polynomials, and, there- 
fore, we give less details. Again, three different cases arise, depending on 
/i = 0, —1, —2, ...,[i = l,2,3,..., and [i otherwise, real or complex. In the 
first and third case we use the Cauchy integral 



E£(z) = ^- / , %r , (3.1) 



where C is a circle around the origin, with radius less than ir. This follows 
from (1.2). 

3.1 Asymptotic form when [i = 0, —1, —2, . . . 

We proceed as in §2.1 and write [i = —m (m = 0,1,2,...). We expand 
En(z) in terms of a finite sum. We have 

££ m (*)=2- m £(T )(z + rT. (3.2) 

For any given m G N and complex z only the term or terms with the largest 
values of \z + r\ will give a large contribution to the sum in (3.2), the other 
terms being exponentially small in comparison. We conclude, that (3.2) gives 
the asymptotic form when n — > oo, when m and z are fixed. In particular, 
when z > 0, the term with index r = m is maximal, and we have 



E~ m {z) = 2~ m (z + m) 



l + O 



z + m — 1 
z + m 



(3.3) 



The error term can also be estimated by C(exp(— n/(z + m))), which is 
indeed exponentially small compared with unity. 
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For general complex z = x + iy and x > —m/2 the term with index 
r = m again is maximal and the same estimate as in (3.3) is valid. When 
x = —m/2 the terms with r = and r = m give the maximal contributions, 
and we have 



E~ m {z)~2- m {-\ m + iy) n + {\ m + iy) n 



(3.4) 



When x < —m/2 the term with index r = is maximal, and we have 



E~ m {z) = 2~"V 



l + O 



z+l 



(3.5) 



As explained at the end of §2.1 these estimates can also be derived by 
using the saddle point method. 

3.2 Asymptotic form when /x = 1, 2, 3, . . . 

We write fj, = m. For m = 1 we have 

e (2k+l)wiz 



E 1 n (z) = E n {z) = 2n\ Y, 



((2k + l)m) n + 1 

E° sin((2/c + Vjirz — \ 
. ±_ 
((2^ + 1)^)"+! 



irn) 



(3.6) 



k=0 



where z G (0, 1) if n = and z € [0, 1] if n > 0. This expansion follows from 
(3.1) as in §2.2. 

In the second series in (3.6) only the term with k = is relevant for the 
asymptotic behavior, and we obtain for fixed complex z 



^) = ^5r[sin (vrz - \ m ) + O (3" B ) 



n — > oo. 



(3.7) 



For general fixed real or complex z the series in (3.6) can be viewed as 
asymptotic expansion for large n, as easily follows from the ratio test. 

For general [/, = m = 1, 2, 3, . . . a similar expansion as in (3.6) can be 
given. In that case the poles of the integrand in (3.1) are of higher order. 
We can write 



3 (2fc+l)7ri2 



k {(2k + l)7ri) n+1 
where e\(n,z) = l,Vfe. 



(3.8) 
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To obtain e™(n, z) for m > 1 we compute the residues of the poles at 
(2k + l)m of order m of the integrand in (3.1). For this we compute the 
coefficient d m -\ in the expansion 



(w - (2k + l)m) m e z 



(e w + l) m w 



n+l 



Y2<I r (w - (2k + l)Tri) r 



r=0 



We substitute w = s + (2k + l)vri and write the expansion as 



( — \^ m e z(2k+l)ni 



(e s - l) m (s + (2k + l)7ri) n+1 



r=0 



(3.9) 



(3.10) 



We use (1.1) and conclude that 



dm-l 



(-l)m e z(2k+l)wi g 



-n — 1 

((2ife + l)7n) n+1 ^ i/! lm-1 - v 



:\v+l—m 



It follows that 



ef(n, z) = (-ir-2- 1 £ ^ ( 



((2k + 1)tt*) 

(3.11) 

\((2k + l)m) v+l - m , 



which we write in the form 

2 m-l 



e£>,z) 



m— 1 



((2/c + 1)ttI 



,m— 1 



n + m — 1 
m — 1 



E 



!/=0 



m — l\ (n + m — z/ — 1) 



v ) (n + m — 1)! 



■(-(2fe + l)7ri) I/ . 



(3.12) 



(3.13) 



For large n the main term occurs for u = 0, giving 

2 



•>m— l^m— 1 



r [l + 0(n^)], (3.14) 



(m- l)!((2fc + l)vri) r 

and in (3.8) the terms with k = 0,-1 give the main terms, and we obtain 
for fixed m and complex z (cf. (3.7)) 

2 m+1 n\ (n + m - V 



7T 



n+m 



m — 1 



TO — 1 



E 



y=0 



"' " ^ " l;! ,"sinr.cO(3 ») 



/ (n + m — 1)! 



(3.15) 



as n — > oo, where r = (z — — ^(m — 1) — ^ iy ) 7r - 
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3.3 Asymptotic form for general complex \x 

The analysis is as in §2.3. We use a contour for the integral (3.1) as in 
Figure 1, now with loops around the branch points ±7u, and with radius of 
the large circle smaller than 3n. We denote the integrals around the loops 
by I±. After the substitution w = 7riexp(s) we obtain for the upper loop 



27ri (7ri) n +M J c 



where 



h(s) = e z 



iris 
e u - 1 



h(s)s-^e- ns ds, 



u = iri (e s — 1) . 



(3.16) 



(3.17) 



We expand h(s) = Xlfc^o ^kS k and interchange summation and integration 
in (3.16). By using (2.33) and (2.34) we obtain the result 



Tr n +v r(/i) 



COS 



(1 - Vhh ( [ ] 



sin 



X Yl 



fc=o k=0 
as n — > oo, with 2 and /U fixed complex numbers (fi ^ Z), where 
X = C - I^tt, C = {z ~ 5^)tt- 



(3.18) 
(3.19) 



The first few coefficients /i^ , are 



/i (r) - 1 
n Q — 1, 



(r) 



^(3(1 - 27r 2 )/i 2 + (13vr 2 - 12Cvr - l)n - 12( 2 
^z{-^ + (1 - vr 2 )^ 2 + 2(^ 2 + 6C 2 )/i - 24C 2 ; 



■>) 
l 2 

dr) 

4 i) = iC(3^ 2 + (vr 2 -7) / x-4C 2 + 4). 
The first-order approximation reads 



cos7r(z — \n — |n) + 0(l/n) 



(0 _ 



(0 



0. 

c, 

(3.20) 
00. (3.21) 
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3.3.1 An alternative expansion 

We repeat the steps of §2.3.1. We write 

'vfl + ir 2 Yf 1 



g(w) 



2ir J V e w + 1 



and expand 



We have 



9(w) =^2("fk + w5 k ) (w 2 + 7r 2 ) fc . 
fc=o 



g(m)+g(-m) 

7o = g = cos ^ 

g(m) - g(-iri) 1 . 

Oo = ^— = -sm(. 

2m tt 

where ( = (z — \h)tt- 

The next coefficients follow from writing go(w) = g{w) and 

g j+1 (w) = SjM-^ + wSj) 



oo 



W 2 + TT 



2 



£ (7k +«><**) K + t 2 ) J \ J = 0,l,2,. 
fc=j+i 



This gives 

= — 



Am' , 

^■(7ri)+^(-7ri)-2(5 J - 
= 4t7 2 ' 



and 



71 = -^[/UCOsC + Tr^sinC], 

5 1 = ^■I 7r? 7 C0S C+(2-M)sinC], 

72 = ^ i [(-9^-37r 2 77 2 + 7rV + 3Ai 2 )cosC + 67r?7(^-l)sinC], 

5 2 = [67rr?(3 - /i) cos C + (36 - 21 fi + 3^ 2 + tt 2 ^ - 3tt V) sin 

96^° 

where rj = \x — 2z. 

Substituting the expansion in (3.23) into (3.1) we obtain 



oo 
k=0 



17 



where 



k-n dw 



w 



n+l 



We have * 



(2n+l) 



and 



V|/(2™) _ ^2k-2ji-2n 



V n 7 1 J n! ' 



(3.29) 



(3.30) 



Hence, 



<(z) = (2^(2n)!^ 7fc ^ 2 " ) , 



fc=0 



< +1 (z) = (2^(2n + l)!^^ 2n) . 

fc=0 



(3.31) 



These convergent expansions have an asymptotic character for large n. This 
follows from 



(2n) 



fc+1 = ,,2 



fi — k 



9 



(2n) 



/i — A; — 1 + n 



(n x ) , n — > oo. 



(3.32) 



Comparing the first term approximations given in (3.21) and those from 
(3.31) we obtain from (3.21) 



ii 



and from (3.31) 

E$n(') = (- 1 ) 



(2n)!2 2 ^ T(n + ^) 



cosvr(z- kt) + ..., (3.33) 



cos ir(z — \n) + . . . . (3.34) 



and we see that the first approximations give the same asymptotic estimates. 

Integer values of fi The expansions in (3.31) do not vanish when (i is 
a negative integer, as the expansion in (3.18) does. We have when \i = m 
(integer) 



(2n) 



( w 2k-2m-2nfk fc > n + ^ 

0, k < n + m. 

Hence, the summation in (3.31) starts with k = n + m. 



(3.35) 
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When we expand 

oo / 2 I 2 \ k 



fc=0 

we obtain the expansions 



k=0 

oo 

< +1 (z)~(2^(2n + l)!^4,^ 2n) , 

fc=0 



(3.37) 



where 7& and 8^ can be obtained from a similar scheme as in (3.26). The 
functions ^jjj 2 "^ are given by 

,( 2n ) _ „-2/i-2n A ~~ M \ _ /_-|\n+fc_-2/i-2n_(^ ~~ ^)n+fc 



*r = -- 2/wn ( v ; + y = (-ir**-^ ^ (n+ y - (3.38) 

When /j, = m (integer) these functions vanish if k — m = 0,1, 2,..., which 
is more useful that in the earlier choice (3.23). For example, when m = 1 
the expansions in (3.37) have just one term (k = 0). In this case ^[j 2 ™'* = 
(— l) n /ir 2n+2 and 70 = 70, <5o = #0 (see (3.24)). These approximations 
correspond exactly to the first term in the second expansion in (3.6). 
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